High Energy Positrons and the WMAP Haze from Exciting Dark Matter 
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We consider the signals of positrons and electrons from "exciting" dark matter (XDM) annihi- 
lation. Because of the light (m^ ^ 1 GeV) force carrier 4> into which the dark matter states can 
annihilate, the electrons and positrons are generally very boosted, yielding a hard spectrum, in 
addition to the low energy positrons needed for INTEGRAL observations of the galactic center. We 
consider the relevance of this scenario for HEAT, PAMELA and the WMAP "haze," focusing on 
light (m^ ^ 2m7r) 4> bosons, and find that significant signals can be found for all three, although 
significant signals generally require high dark matter densities. We find that measurements of the 
positron fraction are generally insensitive to the halo model, but do suffer significant astrophysical 
uncertainties. We discuss the implications for upcoming PAMELA results. 



I. INTRODUCTION 



00 

^vq ■ An overwhelming amount of evidence has estabhshed cold dark matter (CDM) as the standard paradigm for the 
missing matter of the universe. Beginning with early observations of velocity dispersions of galaxies in clusters [l|, and 
later measurements of galactic rotation curves 0,[3|, CDM has been supported by numerous additional observations. 
These include strong lensing of background galaxies 0] , x-ray emission from galaxy clusters ^5|| , the combination of 
I . CMB and type la supernovae data measurements of the distributions of galaxies 0, [1], as well as the highly 
O ' remarkable recent study of the bullet cluster [9i|. 

^ The nature of this dark matter remains an open question. Measurements of the CMB, the need for early structure 

growth, and the success of big-bang nucleosynthesis rule out baryonic matter as being the dark matter, necessitating 
a new particle beyond the standard model. One of the most appealing of these is the thermal WIMP, whose relic 
abundance has a simple relation with the annihilation cross section. Specifically, for thermal s-wave freezeout, one 
^ ; has [3 

■ ni /2.5X 10-26 cm3s-l\, , ,^ , 

CN . nh' = 0.1 X — /(to) (1) 

<N ; 

where /(m) is a logarithmically varying function of mass, with /(500 GeV) — 1. Because {aann = 2.5 x 
. IQ-^^ cm'^s"^ is a cross section appropriate for a weak scale particle, a scale already suggested by the hierarchy 
' problem, we have strong motivation to consider a WIMP with a mass in the range of a few hundred GeV. 

o ■ 



A. Indirect Detection of Dark Matter 

Indirect detection involves observing either the annihilation products of dark matter or the decay products of 
unstable states. These products include positrons, photons, anti-protons, neutrinos, and anti-nuclei. Searches are 
often made for the anti-particles, because their astrophysical backgrounds tend to be smaller than those of particles. 
Though the DM annihilation cross-sections are generally small, if the density of dark matter particles is large, as 
is believed to be the case in the center of the Galaxy, then the annihilation rate can be large enough to produce 
observable effects. Currently experiments are underway to detect neutrinos (AMANDA JlTl, IceCube from 
annihilation of dark matter captured in the Earth and Sun, and to detect anti-particles (PAMELA [l^) and photons 
(GLAST dl, HESS [l5], INTEGRAL 0) from dark matter annihilation in the Galactic halo. 

In many theories of SUSY dark matter, the dark matter particles are Majorana fermions that annihilate through 
processes such as XX — ^ If- Annihilations to heavy fermions are favored over those to the light fermions, because 
the cross-section is proportional to TOj. Therefore, the and e~ spectra are generally soft. This can make the DM 
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signal difficult to detect, as the Galactic background spectra for electrons and positrons are largest at low energies 
and are currently not precisely understood. For this reason, a theory of dark matter that produces a hard spectrum 
is perhaps the best chance for an observable indirect detection signal from electron and positrons. 

Recently, a number of different experiments have made measurements which are consistent with excesses in positrons 
and/or electrons that may be indicative of a dark matter signal. Two of these (HEAT and the WMAP haze) are 
indicative of high energy particles, whereas the results of INTEGRAL/SPI are indicative of low energy particles. 

In this paper, we will explore the possibility that a single source can produce all of these signals, specifically the 
"exciting dark matter" (XDM) proposal. We will focus on the cases most favorable for achieving high energy positron 
signals, in which the annihilation products of the WIMP are lighter than about 250 MeV, and decay either into ii'^/J,~ 
or e~^e~ pairs. In the next section, we will review the evidence for excess electronic activity which may arise from 
dark matter. We then will describe the XDM scenario, and how it produces large quantities of high energy electrons 
and positrons. We describe the calculation of observable positron fraction and synchrotron signals in section [TVl In 
section |V] we present our results. We find the XDM can easily explain both the HEAT excess as well as the WMAP 
haze, although significant signals tend to require large densities of dark matter. We find a wide range of spectra 
that may be found at PAMELA, including very hard positron spectra arising from the highly boosted annihilation 
products of XDM. Finally, in section IVTl we conclude. 

II. SIGNALS FROM POSITRONS AND ELECTRONS 

Positrons and electrons are the products of dark matter annihilation in numerous theories of WIMP dark matter, 
so their careful study could prove fruitful in the quest to understand the nature of dark matter. The existence of 
electrons and positrons can be observed directly through a measurement of their particle fluxes, or indirectly through 
the measurement of their associated radiation, including 7-rays resulting from their annihilation and synchrotron 
radiation. 



1. Direct Measurement of e'^ and e Spectra 

Direct measurements of e"*" and e~ spectra have been done in several balloon-borne cosmic ray experiments, most 
recently in HEAT2000 and CAPRICE98 PAMELA, a satellite experiment, is currently taking data. (Note: Many 
experiments state their results for the ratio e+/(e~'" + e~), the positron fraction, to eliminate some systematic errors 
in the measurements of the individual fluxes.) Cosmic ray electrons and positrons are of two varieties: "primary" 
particles created by astrophysical sources, and "secondary particles" produced through the interactions of primaries 
with gas in the interstellar medium (ISM). Above 100 MeV the primary positron flux due to standard astrophysical 
sources is thought to be negligible [18], while the secondary flux behaves like an inverse power law. Therefore, the 
positron background signal is small at large energies, which makes the positron signal a promising candidate for 
observing new physics. 

For instance, the balloon-borne HEAT (High-Energy Antimatter Telescope) experiment flew twice in the mid-1990's, 
and measured the individual and combined energy spectra of electrons and positrons [l9| . From these measurements 
the positron fraction was determined for the energy range from 1 to 50 GeV. The HEAT results suggest that for 
energies larger than 10 GeV an excess in the positron flux above that expected from pure secondary production 
cannot be ruled out [T^ . There may exist sources that give rise to primary positrons with energies above 10 GeV. 
Additionally, there appears to be some structure in the positron fraction above 7 GeV that a pure secondary spectrum 
cannot explain [20| . 

An important experiment which will clarify the results of HEAT is PAMELA. The Payload for Antimatter Matter 
Exploration and Light- nuclei Astrophysics experiment is a satellite-borne apparatus launched in 2006 that is currently 
measuring the energy spectra of many cosmic ray particle species, including electrons, positrons, protons, anti-protons, 
and light nuclei. One of the stated primary objectives of the experiment is to search for evidence of annihilations 
of dark matter particles by measuring the electron and positron energy spectra in the energy range from 50 MeV to 
270 GeV [13]. PAMELA will extend the measurement of the positron fraction out to an energy almost an order of 
magnitude larger than that achieved by the HEAT experiment. Additionally, the PAMELA results will be based on 
very high statistics, approximately 10^ e+ per year [ij]. Initial high energy results from the PAMELA experiment 
are expected soon. 
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2. High Energy Positron and Electron Signals in Synchrotron Radiation 

In addition to directly detecting the cosmic ray positrons, one can detect 7-rays coming from electron-positron 
annihilation and synchrotron radiation coming from electrons and positrons propagating through the galactic magnetic 
field. An excess of either of these radiation signals could indicate the existence of dark matter annihilations. The 
former is generally a measure of stopped or low energy positrons, while the latter is a measure of higher energy 
(multi-GeV) positrons and electrons. 

Let us begin by considering the possibility of synchrotron radiation. The Wilkinson Microwave Anisotropy Probe 
(WMAP), launched in 2001, measured the microwave emission at five frequencies, 22.5, 32.7, 40.6, 60.7, and 93.1 
GHz, over the full sky. Multiple frequency measurements were necessary to separate Galactic foreground signals, each 
with a unique spectrum and spatial distribution, from the Cosmic Microwave Background (CMB) signal. 

Using the one-year data from WMAP released in 2003, Finkbeiner identified the microwave emission due to the three 
well-understood Galactic signals, thermal radiation from dust, free-free (bremsstrahlung) radiation, and synchrotron 
radiation from high-energy electrons (accelerated by supernovae shocks) spiraling in the Galactic magnetic field, and 
additional emission due to spinning dust [2l|. His analysis indicated that there is excess emission, not correlated 
with any of the known Galactic foregrounds, that is distributed with approximate radial symmetry within ~ 20° 
of the Galactic center and that decreases rapidly with projected distance from the Galactic center. He called this 
excess microwave emission the "haze" and argued that it is consistent with synchrotron emission from high-energy 
electrons and positrons created in the center of the Galaxy [l^] . One explanation for the presence of these electrons 
and positrons in the inner Galaxy is that they are produced by annihilating dark matter [l^, [l^l ■ 

3. INTEGRAL/SPI 

The strongest gamma-ray line signal from our galaxy comes from electron-positron annihilation [24]. The SPI 
imaging spectrometer aboard ESA's INTernational Gamma-Ray Astrophysics Laboratory (INTEGRAL) gamma-ray 
satellite observatory has measured the 511 keV line coming from direct annihilation of e^e~ pairs and decay of 
para-positronium, as well as the continuum spectrum from the 3-photon decay of ortho-positronium. Both signals 
are strongest within a few degrees of the Galactic center, indicating that this region has the largest concentration 
of electron-positron annihilation [23|. The 511 keV emission from the bulge has ellipsoidal symmetry about the 
Galactic center with an angular extension of 6.5°tJ;g (FWHM) in longitude and 5.1°;p;8 (FWHM) in latitude [H]. 
The flux from the bulge at 511 keV is (7.04 ± 0.32) x lO^'' ph cm^^ s~^, while the emission from the disk is 
(1.41 ± 0.17) X 10-3 ph cm-2 s-i [H. 

It is presently a great challenge for conventional astrophysical sources, such as cosmic ray interactions with the 
ISM, neutron stars, black holes, supernovae, low mass x-ray binaries, and pulsars, to explain the INTEGRAL signal, 
although, with myriad uncertainties, they may prove to be responsible. Recently, it was noted that the disk component 
of positron emission, in particular, was likely due in large part to LMXBs [25| . although it is still uncertain whether 
such objects can provide a significant piece of the bulge emission. 

A very appealing possibility is that the excess arises in some fashion from dark matter. However, such explanations 
are challenging, because the INTEGRAL signal is inconsistent with injected positrons at energies much higher than 
a few MeV (2^, [l^, HI] . As such, we cannot simply identify these positrons with annihilation products of weak-scale 
dark matter. 



III. EXCITING DARK MATTER (XDM) 

Theories of "light" dark matter, for example scalar dark matter [29l[30l[3]|. with masses in the MeV range well below 
the weak scale, as well as decaying dark matter scenarios [3l,|33|, have been proposed to explain the INTEGRAL signal. 
However, these theories are ad hoc to some degree from a particle physics perspective. An alternative explanation 
for the INTEGRAL signal with a weak scale dark matter particle arises in the scenario of "exciting" dark matter 
(XDM) [33|. XDM proposes that the electron-positron pairs needed to explain the INTEGRAL signal are created in 
the decay of an excited state of WIMP dark matter. The theory contains a weak scale dark matter particle x with 
an excited state x* that has an energy at least 2me above that of the ground state. The coupling of the pseudo-Dirac 
dark matter particle to standard model particles is through a light scalar particle 0, ~ .1 — 1 GeV, which couples 
to the Higgs. The decay of the excited state into the ground state can result in the emission of an e+e" pair. 

The energy needed for the excitation process comes from the kinetic energy of the WIMPs; a 500 GeV particle 
with a speed of 500 km/s has a kinetic energy greater than me. XDM can account for the creation of 3 x 10''^ e+e" 
pairs per second in the Galactic bulge, enough to explain the INTEGRAL 7-ray signal. Although the model requires 
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large scattering cross sections, and cuspy halos, it is an intriguing suggestion that allows one to generate low energy 
positrons from weak-scale dark matter particles. Beyond particles produced from excitation and subsequent decay 
of WIMPs, additional signals of XDM arise through annihilation of x iiito 0, which proceeds through the processes 
shown in Fig. [TJ 




FIG. 1: Diagrams contributing to annihilation of x- 

Since present-day dark matter is moving very slowly (/? ~ 10"'^), only s-wave processes are relevant for indirect 
detection in the halo (at least for thermal dark matter). Since we are motivated by possible existing signals, it is 
important to consider whether XDM can produce models with significant s-wave annihilations. This is not a trivial 
point. Although XDM functions as a Dirac fermion at freezeout, where the splitting between x ^^nd x* is small, 
in the present universe, by assumption, only x is present. This is a Majorana fermion, and thus to annihilate via 
s-wave, must be put into a state with CP — -1 (see the discussion in [1^). As a consequence, Majorana XDM cannot 
annihilate into two identical real (/) particles. However, relatively minor modifications can allow this annihilation to 
proceed via s-wave. In particular, we can allow 0(1) CP violation in the dark matter sector, or we can promote (j) 
to a complex scalar. In models of XDM where the dark matter is a Dirac fermion with an excited state (in analogy 
with hadrons), no such modifications are necessary. In any event, the modifications to the model in [34] needed to 
generate s-wave annihilations are small, so the resulting differences in the signals in question are small. 

Since ^ to^, the (j) particles come out boosted by a factor of 7 = The (/)S can then decay into SM particles 
through their mixing with the Higgs. For < 2me, 4> can only decay into two photons. For 2me < < 2m^, the 
decay can proceed through two channels, (j) — ^ 77 f-nd 4' — * e~e"'". The decay of 4> to two 7s is due to one loop 
spin-0, spin-I/2 and spin-1 contributions. These almost cancel giving a branching ratio (B.R.) for — > 77 of at most 
14% in that energy region j36l |. In the energy range 2r7i^ < < 2-01^^0, the channel — > M /^^ becomes relevant, 
and since the branching ratio of 4> to // goes as m?j, the decay (p — > /^"M^ completely dominates over (j) — > e^e"*". 
In this energy range the decay of (p into 2 photons is negligible [s^ . Above the threshold energy 2to^o , the channel to 
neutral pions is opened. For (f) masses above 2m^±, it has been shown that the B.R.s for the processes (j) — * 7r~7r+ 
and (j) — > 27r° cannot be ignored [13, IH, [3|. For a 1 GeV (j) the B.R. for — > ^-^+ is w 20% ^ while for the 
decays to tts it is « 40%. For > 1 GeV decay channels to heavier mesons appear. The work presented in this 
paper investigates the mass range 2me < < 2mTr. We leave other mass ranges for future work. 

We considered two distinct cases for the mass of m^: = 0.1 GeV, for which the B.R. for — > e^e+ is = 90%, 
and = 0.25 GeV, for which the B.R. for decay to and /i^ is = 100%. In both scenarios, the resulting spectra for 
the electrons and positrons are much harder than typical e+e^ spectra coming from weak scale WIMP annihilation. 
The positron fraction has a characteristic bump that starts around 6 GeV. (The precise location of the peak depends 
on the masses of x ^-nd 4>.) There is evidence for additional structure around this energy in the HEAT measurements 
of the positron fraction; XDM provides an explanation for this. The presence of additional high energy e^e~ pairs 
from XDM annihilation in the center of the Galaxy would give rise to additional synchrotron radiation as they spiral 
through the Galactic magnetic field. In this way XDM provides an explanation for the excess microwave emission, 
the "haze" [Hi, observed by WMAP. 

IV. CALCULATIONS 

We have calculated the local positron fraction e+/(e+ + e~) up to energies of 1 TeV and the local synchrotron 
radiation due to dark matter for the haze frequencies, 22.5, 32.7, 40.6, 60.7, and 93.1 GHz, using GALPROP ^4^ 
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version 50p j72| . The ratio calculation requires knowledge of the local primary and secondary positron and electron 
spectra, as well as the contributions to both local spectra from dark matter particle annihilations. The contribution to 
the positron fraction from XDM annihilation is significant only for energies above ^^5 GeV, so our starting assumption 
is that the background fraction should fit well the HEAT data at energies below '^5 GeV. We calculated our background 
ratio in agreement with the first five HEAT data points and required the total positron fraction to fit well all nine 
HEAT data points. In light of prelimary PAMELA data on positrons below 10 GeV, we also fit to those, as well. 
We predict the results of the PAMELA experiment for the positron fraction in the XDM scenario. Additionally, we 
compare our calculation of the XDM synchrotron radiation to the "haze" . 



A. The Injection Spectra for e e Arising from DM Annihilation 



The annihilation process for XDM is through the intermediate 4> particle. The (j) mixes with the Higgs allowing it 
to decay into the final state fermions. We consider two scenarios. The first is the "direct" decay channel, in which 
the annihilation proceeds directly into e~ pairs, which occurs for ra^ < 2m^. For 2m^ ^ ^ 2m^o we have the 
"muon" decay channel, where 4> decays to muons, which in turn produce electrons and positrons. Higher mass (f) can 
decay into pions, kaons, and other hadrons. We defer the analysis of these scenarios to a future paper. For the direct 



channel, the injection spectrum is flat (see Fig. |2(a) ), while for the muon channel it is somewhat softer (see Fig. 2(b) I 
We give the details and formulae in Appendix [Al 
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FIG. 2: Injection spectra of positrons. 

It is important to note that these light particles produce few 7-rays directly, and are kinematically incapable of 
producing anti-protons, which are a strong constraint on these models. The absence of indications of anti-proton 
excesses may be indicative of such a light mediator producing a positron signal. 

The GALPROP code of Moskalenko and Strong calculates cosmic ray propagation through the Galaxy by solving 
the propagation equation numerically on a grid [40| • The spatial dynamics include diffusion resulting from cosmic rays 



scattering on MHD waves, and convection on Galactic winds 41|. In momentum space, diffusive re-acceleration, the 
result of stochastic acceleration due to scattering on MHD waves, and energy loss from ionization, bremsstrahlung, 
inverse Compton scattering and synchrotron radiation are included [iO, |4l[ . Using the most recent measurements of 
the source abundances for the primary species (nuclei, electrons and 7-rays), GALPROP propagates the primaries 
through the Galaxy, then iteratively computes the resulting spallation source functions for all species, and propagates 
the equivalent full, primary + secondary, source function for each species until a converging result is obtained [i^ ]. 
GALPROP assumes free escape of the particles as the spatial boundary condition [isj. For our calculations, we used 
a 2-dimensional spatial grid (and a one dimensional energy grid), and assumed cylindrical symmetry as well as mirror 
symmetry with respect to the Galactic plane. 

Through the galdef file, the parameter input file, the code allows for considerable freedom in the choice of many 
astrophysical parameters, including the primary electron injection spectrum, the diffusion coefficient, and the strength 



of re-acceleration. We used the "conventional" model of GR production and propagation as used by Abdo et al. [44 1 



as the starting point for our choice of parameters. The conventional model calculates the local proton and electron 
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TABLE I: Parameters of the primary electron spectrum for fits to HEAT 



VA 

(km/s) 


Injection Index 
Below Break Rigidity 


Break Rigidity 
(MV) 


Electron Flux Norm. 
(cm-^sr-^s-^MeV^^) 





1.60 


4.0 X 10^ 


0.2488 X 10"" 


20 


2.10 


6.0 X 10^ 


0.2612 X 10"'' 


35 


2.15 


4.0 X 10^ 


0.2887 X 10"" 



spectra in agreement with the locally measured values f45'| . For a given choice of Alfven velocity, we varied only four 
physical parameters to fit the low energy HEAT data, specifically, the energy of the break in the primary electron 
spectrum (Break Rigidity), the injection index below and above the break, and the overall flux normalization (see 
Table HI). We also vary one calculational parameter, the energy grid spacing, in our calculations of the positron 
fraction. 

B. Dark Matter Density Functions (Halo Profiles) 

In order to understand the dependence of the signals upon halo models, we performed our calculations using three 
different halo profiles: 



1 



N.F.W. Profile 



Pir) = Po 



J.2 _|_ J.1 



Isothermal Profile 140 



2 r" - 

p{r) = Po exp[ ( ; 

a r 



Merritt Profile m 



Here Rq = 8.5 kpc is the solar distance from the Galactic center, = + is the spherical radial coordinate, 
.13 <a< .22 [t^ is the Merritt parameter which defines the cuspiness of the profile, and r_2 = 25 kpc is the radius 
at which the logarithmic slope of the Merritt profile is -2. Tc is the core radius and po is the local value of the dark 
matter mass density. The values of these parameters for the three profiles are listed in Table HT] See Fig. [3] for a 
comparison of these three DM density profiles. We limited ourselves to these three halo models, which arise under 
various assumptions and circumstances, and did not consider theoretically unmotivated profiles, such as a fiat halo 
profile. The Merritt profile, being the cuspiest, is generally the best for generating the INTEGRAL signal from XDM. 
The isothermal profile provides a theoretical simplification, which is often employed for analyzing, for example, direct 
detection limits. 

TABLE 11: Parameters of the dark matter density profiles 



Model 



(kpc) 



N.F.W. 

Isothermal 
Merritt 



20.0 
2.8 
25.0 



C. The Electron and Positron Backgrounds 



The background positron fraction is given by 



secondary 



^secondary ~^ ^primary ~^ ^secondary 
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FIG. 3: Dark Matter density profiles. 



The primary positron spectrum is taken to be zero, as it is negligible in comparison with the secondary spectrum. 
The sources of the primary e~ spectrum are believed to include supernovae remnants and pulsars, though their exact 
nature is unclear. Additionally, there is uncertainty in the primary e~ spectrum itself, particularly at energies below 
10 GeV, where the effect of solar modulation is to attenuate the spectrum at lower energies. However, it is generally 
agreed that for the energy range of interest to us, 1 GeV to 1 TeV, the spectrum is well described by a power law 
E~°' with an index a that increases at higher energies. The GALPROP code provides the freedom to modify the 
primary electron injection spectrum. Inputs into GALPROP include constant values of a in three energy regions, the 
two energies defining these regions, the normalization of the electron flux, and the energy at which this normalization 
holds. 

Secondary positrons and electrons are created in the interactions of primary CR protons and Helium as they 
propagate through the Galaxy, p-p and p-He collisions result in charged pions and kaons, which give rise to electrons 
and positrons through the decays K"^ — > tt^ + 7r°, ^ fj,"^ + Vf^, tt^ ^ + f^, and ^ + iTJ^ + Ve- The 
secondary positron and electron spectra are equal to one another by CP invariance of weak decays [i^ . 

We used a standard least-squares fit of the calculated background ratio to the first five HEAT data points. The 
backgrounds are shown in Fig. [H It is clear from the figure that the backgrounds are very sensitive to the choice of 
parameters for the primary electron spectrum. In spite of the differences in the backgrounds, as we show later, there 
are common features in the total positron fraction. As a consistency check on our backgrounds, we calculated the 
resulting local H and He fluxes and compared them to the measured fluxes at the top of the atmosphere [49, 50, 51,]. 
We had good agreement for both spectra for energies above a few hundred MeV. It is important to note that differences 
in the low energy background carry over to high energies, and thus lead to differences in the positron fraction at high 
energies. 



D. The Positron Fraction 



We calculated the total positron fraction ^] for each of the three backgrounds, va — 0, 20, and 35 km/sec, 
varying the following parameters as described: five dark matter halo profiles, N.F.W., Isothermal, Merritt (a = 0.13), 
Merritt(a = 0.17), and Merritt(a = 0.22); five DM masses, m^= 50, 100, 250, 500, and 800 GeV; and two decay 
channels, direct and muon, corresponding to = .1 and .25 GeV, respectively. We fit the total fraction to all nine 
points of the HEAT data using a least squares fitting routine with the annihilation cross-section as the fit parameter. 
As stated previously, a value of the thermally averaged annihilation cross-section {aann \v\) that gives the correct DM 
relic abundance is roughly 2.5 x 10~^^ cm^s~^. However, because the physical signal, in particular for high energy 
positrons, is dominated by the local density, we quote our results in terms of the physically relevant parameter (p^crw). 

Inhomogeneities in the dark matter density can lead to additional "boost" factors, b = (p^) / (p)'^ ■ Based on 
present halo simulations, a reasonable value for a boost factor for gam ma ray production in the halo is 5 ^ 3, although 
extrapolations to smaller scales maybe allow boosts as large as 13 [52|. Here, we are interested in two separate signals, 
both from the inner region of the halo, and it is not yet clear what a reasonable boost factor is. However, it is more 
likely that clumps are tidally destroyed in the region closer to the galactic center, thus we expect, for reasonable 
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Positron-Electron Ratio Bacl<grounds 



4 = km/s, X . 6.61 ■ 

- 20 km/s, . 3.97 . 

= 35 km/s, ^ - 1 .30 ■ 

HEAT Data i 



Energy (GeV) 

FIG. 4: Positron fraction backgrounds for va = 0, 11,4 — 20, and va = 35 km/s. The quoted is for the fit to the five lowest 
energy points. 

scenarios, that a larger boost factor would arise for the positrons produced away from the GC, i.e. those relevant 
for HEAT and PAMELA, than for positrons produced near the GC, i.e. those relevant for the haze. That said, 
uncertainties in the galactic magnetic field and halo profile can also play an important role in the haze signal, so the 
boost factor that we report here is still uncertain. 

E. DM Synchrotron Radiation 

The calculation of the synchrotron radiation due to e~ pairs arising from dark matter annihilation is unencum- 
bered by the complications of the background fluxes, since the dark matter contribution is completely independent 
of the contributions from standard astrophysical sources. We fit the calculated intensity for 22.5 GHz synchrotron 
radiation to the 29 haze data points extending over 6° - 34° from the Galactic Center using the annihilation cross- 
section and a constant offset as fit parameters. (The inner 5° are masked due to the bright dust emission.) Because 
we only expect dark matter to dominate roughly the inner 15°, and because the overall normalization of the haze is 
somewhat uncertain, we fit the existing data to a dark matter component, plus a constant offset. This, then, prevents 
the 15° — 35° region from dominating the fit. 

F. Uncertainties 

A given model generally makes a very robust prediction for the primary positron injection spectra once all kinemat- 
ical elements are specified (specifically, dark matter masses and the masses of the annihilation products). However, 
when converting these injection spectra into observed rates and spectra at an experiment such as PAMELA, there 
are many uncertainties that must be considered. 

1. DM Halo Uncertainties 

The amplitude of the signal is first and foremost determined by the cross section and the density of the dark matter 
pa, entering in the form Pq (a \v\). While a thermal WIMP will generally have a cross section 2 — 3 x 10~^^cm'^s~^, 
non-thermal WIMPs can have significantly higher cross sections giving higher rates at HEAT and PAMELA and for 
the WMAP haze [Ulllllil- 

The local dark-matter-produced positron signal at the highest energies is simply determined by the local density 
of dark matter. This is because high energy positrons propagate a very short distance before losing much of their 
energy. A 300 GeV positron, for instance, typically drops to x 300 GeV in a distance of ~ 300 pc. Consequently, 
the signal for HEAT and PAMELA is relatively insensitive to the broad shape of the halo, but is very sensitive to the 
local value of po- 
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Conventionally, po = 0.3 GeVcm"^ is adopted for many analyses. While providing a useful baseline, it is important 
to recognize that there is significant uncertainty in this value. The standard reference for a determination of pq 
is [HBl. There, the expected range is given as 0.2 GeVcm"'^ ^ GeVcm^'^, while noting simultaneously that 

additional flattening of the halo could increase this by a factor of two. (57|] considers halo models with densities up 
to 0.8 GeVcm^"^, while 58] considers local densities of 1.07 GeVcm"^. 

The dumpiness of the halo (commonly referred to as a "boost factor", b =< > / < p >^), is simply the 
observation that local overdensities can increase the annihilation rate. Studies of this for gamma rays have indicated 
that boosts 6 « 3 are quite reasonable, and boosts 6 « 13 might occur, if the boosts found in the simulation are 
extrapolated to the subhalos [s^l, although it is not clear that such boosts survive to interior regions of the Galaxy 
as are relevant for HEAT and PAMELA. Moreover, such overdensities can impact the spectrum if a particularly 
significant overdensity is located near us [ssj . 

Ultimately, it is reasonable to consider values as large as PQ{f7\v\) ~ 3 x 10^^^GeV^cm~'^s~^ without invoking 
boost factors, and possibly significantly higher, p^ {a \v\) ~ 4 x 10~^^GeV^cm^'^s^^, if one considers boost factors as 
well. Non-thermal candidates can invoke yet higher values. We shall find that most regions of parameter space that 
produce significant signals at PAMELA are on the high end of this range. However, in the event of a significant signal 
at PAMELA, it will be unclear whether a local overdensity or boost, a high value of po (perhaps from a flattened 
halo), or large non-thermal cross sections are responsible for it. 

2. Uncertainties in the Propagation Parameters 

Relevance to the Positron Fraction 

We require our results to be consistent with local CR measurements, since it is the local positron fraction we are 
interested in. Local cosmic ray measurements allow for uncertainties in the primary electron and secondary electron 
and positron fluxes, and these in turn give rise to different shapes for the positron fraction. The uncertainties in the 
CR spectra can be associated with uncertainties in the production mechanisms for the spectra, i.e. the shapes of the 
spectra themselves at production, or with uncertainties in the CR propagation parameters. 

The total (primary-f-secondary+dark matter) electron spectrum is constrained up to 100 GeV by measurements 
of the BETS [M], HEAT CAPRICE [6l|, and MASS [62] experiments, among others. Since the secondary and 
dark matter contributions to the spectrum are at least an order of magnitude smaller than the primary contribution 
at all energies, the fit of the total electron spectrum to the data is determined almost exclusively by the shape 
of the primary spectrum. For a given set a propagation parameters, the data determines the following electron 
flux parameters (though it does provide for some freedom in these parameters): the primary electron spectrum 
normalization, power law indices, and power law break rigidity. The values of the injection power law index for 
energies greater than a few GeV that are consistent with local measurements include the values of 2.54 and 2.6, which 
we have used in our background models. Our results deviate from the local electron spectrum at energies below E = 5 
GeV where solar modulation plays an important role. This is because our calculations using GALPROP don't take 
into account the effects of solar modulation. This disagreement is expected to show up as in [4l|. Changes in the 
propagation parameters may give rise to changes in the electron flux. As an example, electrons lose energy relatively 
quickly as they propagate, so the high energy contribution to the local flux must come from electrons produced in 
nearby regions of the Galaxy. Therefore, changes in the propagation parameters that affect propagation on short 
distance scales can affect the flux at high energies. In this way, the uncertainties in the propagation parameters can 
lead to uncertainties in the local particle fluxes. 

Secondary positrons and electrons are created as primary CR protons and Helium propagate through the Galaxy. 
The local proton flux is constrained from a few tenths of a GeV up to ~ 180 GeV by the measurements of BESS [49j . 
IMAX [50], and AMS [5l|. The primary nuclei spectra are described by a power law E~°', with an increasing value of 
a at higher energies. GALPROP allows for defining the constant values of a in two rigidity regions, the value of the 
rigidity that defines those regions (break rigidity), the normalization of the proton differential flux, and the proton 
kinetic energy at which the normalization is being carried out. The proton data allows for some, though not much, 
uncertainty in these proton flux parameters. As a result, there is some corresponding uncertainty in the secondary 
positron spectrum. For energies higher than E ^ 10 GeV, values of a in the range = 2.6 — 2.7 are consistent with 
the locally observed spectra. Lower values, as in [3], where an index value of 2.42 was used, are also within the 
current uncertainties. Higher values of a for energies £; > 10 GeV resuh in a decrease in the flux of the secondary 
electrons and positrons. This in turn allows for a higher (up to a factor of = 2) I ^tecondary * given local DM 
density po and annihilation cross section. Additionally, since positrons and nuclei propagate differently, the choice of 
propagation parameters can affect the positron spectrum for a particular choice of proton spectrum parameters. 

The propagation of cosmic rays is governed by the transport equation, which includes the effects of diffusion and 
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energy losses, among other things. We chose to hold fixed the energy loss mechanisms provided in the GALPROP 
code. However, in addition to varying the parameters associated with re-acceleration, we varied the following diffusion 
parameters: the diffusion zone length 2zmax (the perpendicular distance from the Galactic plane at which free escape 
of cosmic rays is assumed), and the energy-dependent diffusion coefficient, specifically through its normalization, 
break rigidity, and power law index above and below the break. 

Mechanisms for re-acceleration of charged particles within the halo are well motivated, but still poorly understood. 
In order to provide some sense of the effects of this uncertainty, we employed three different values of Alfven velocity, 
a parameter describing the strength of diffusive re-acceleration. (The Alfven velocity relates the spatial diffusion 
coefficient to the momentum-space diffusion coefficient.) In our calculations we take va = km/s corresponding to no 
re- acceleration, and va = 20 km/s and 35 km/s corresponding to the inclusion of re-acceleration in CR propagation. 
These values have been used extensively in the literature [g^ and are within the acceptable range of values ^64, 69| . The 
effect of re-acceleration is to shift the flux of electrons and positrons to higher energies. In our calculations, significant 
contributions to the fluxes occurred only for energies well below 20 GeV, nonetheless, we often see differences in the 
high energy positron fraction when varying va- This arises as follows: our primary and secondary fluxes in the low 
energy range of the HEAT data are highly dependent on the value of the Alfven velocity. In order to fit our three 
backgrounds to the first five HEAT data points, as previously described, we varied the primary electron spectrum 
through the low energy power law index, the break energy of this index, and the normalization. (The values of these 
parameters for the three backgrounds are listed in Table ID) The changes in the low energy background then carry 
over to high energies, and thus lead to variations in the positron fraction at E^20 GeV. 

The GALPROP code provides for an energy-dependent spatial diffusion coefficient D^^ defined as: 

/ R X^" 

D:,^ = PDo^x (2) 

\ ^rigid br J 

where i? is the rigidity of the particles, (3 = v/c, Dq^x is the diffusion coefficient divided by /? at rigidity Dri„id br, and 
the index Dg can take the value Dgi for R < D^igid br and a different different value Dg2 for R > Drigid br j43l |. For a 
given Zmax and re-accleration strength, the diffusion coefficient is determined by fitting the Boron-to-Carbon (B/C) 
ratio to local measurements. In general, an overall decrease in the diffusion coefficient gives an enhancement in the 
flux of low energy electrons and positrons, while having little effect on the proton spectrum. This effect provides for 
additional uncertainty in the secondary positron spectrum; the flts to the proton data can be maintained for many 
variations on the secondary positron spectrum. 

In the models presented here we use Dq^x = 5.8 x lO^^cm^s"^. Increasing the value of Di^^x results in an increase in 
the ratio e^^j^J e^^^^^j^ at the energies of interest to us, 100 GeV, but for values of Dq^x ?i 6 — 7 x IO^^cto^s"^ 
with Drigid br — few Ge V , serious inconsistencies arise between the calculated electron flux and the local electron flux, 
unless we choose a smaller power law index for the primary electron spectrum for energies above a few GeV. 

For energies where the XDM positron contribution becomes comparable to the secondary positron flux, the power 
law index Dg is signiflcant. Values between Dg — 1/3 (corresponding to Kolmogorov turbulence) and Dg — 1/2 
(corresponding to Kraichnan turbulence) are considered reasonable [43l . [6a | . Increasing Dg from 1/3 to 1/2 can increase 
^DM/^tecondary '^^P ^ factor of 2. Howcvcr, to be consistent with the local proton and CR nuclei measurements, the 
primary proton power las index must be decreased in the region of energies higher than E = lOGeV. This cancels the 
gain in e+4^/e+^„„^^^^ by about the same factor. 

The diffusion zone width 2zmax can take on a range of values. The allowed values for propagation models that include 
re-acceleration and for those models with both no re-acceleration and no convection are in the range 4 < Zmax < 12 
kpc [G^I- Increasing the diffusion zone width allows fewer particles to escape from the conflnes of the Galaxy, so 
particles from farther away are able to make their way to our location in the Galaxy. As mentioned earlier, electrons 
and positrons lose energy rapidly as they propagate. The result is that the flux of low energy electrons and positrons 
is enhanced. The positron flux at energies E ^ lOOGeV is affected negligibly, as the highly energetic positrons we 
measure come from a sphere with a much smaller radius than the diffusion zone width. 

Relevance to Synchrotron Radiation (The Haze) 

Since the emissivity of synchrotron radiation at a speciflc frequency is proportional to the magnetic fleld, the most 
dominant parameter determining the shape of the synchrotron spectrum, apart from the DM halo proflle, is the 
magnetic field. A reasonable value for the local magnetic field is i? = 3 with 5 fiG being within the limits. 
For the center of the Galaxy, values within the range of 5 — 20 /iG are acceptable. Little can be said about the 
parametric dependence of the intergalactic magnetic field on z, the vertical distance from the galactic plane, or on 
the radial distance R from the center of the Galaxy at z = 0. We considered, as in |43|, a magnetic field of the form 

B = Bie «<: e , where Bi is the local value of the B-field, and Rc and Zc are the characteristic scales of the 
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B-field, taken to be 10 kpc and 2 kpc, respectively. We took the local value of the magnetic field to be 5 ^G. More 
homogeneous magnetic fields have been used in the literature . 

The synchrotron spectrum also depends on the CR diffusion parameters. A smaller diffusion coefficient suppresses 
the diffusion of e"*" and e~, thus they remain inside the diffusion zone for longer, yielding larger synchrotron fluxes. 
Because the synchrotron radiation at the Haze frequencies gets it's dominant contribution from and e~ with 
energies above a few GeV, the dependence of the diffusion coefficient on the rigidity at high energies can have a small 
effect on the synchrotron flux. Since the Haze is in the region of = 6 — 15 degrees in latitude or r = 0.9 — 2.3 kpc 
from the center of the Galaxy, varying the diffusion zone width above z = 4 kpc doesn't have signiflcant effects on 
the shape of the synchrotron radiation as a function of latitude. 

Additionally, varying the ratio of the energy density in the magnetic field to the energy density in the radiation 
field Umag /Uracil which is equal to the ratio of energy loss of high energy e+e^ through synchrotron radiation to that 
through inverse Compton scattering, has an affect on the synchroton radiation produced [54| . A larger value of Urad 
results in faster depletion of the number of high energy e^e~ pairs, while lower values can yield a stronger haze signal. 

V. RESULTS 

A. Baseline Models 

In order to determine the sensitivity of the signal to various unknown variables, such as particle mass and halo 
model, we consider a particular baseline or "canonical" model about which to vary. We have two canonical models, one 
without and one with re-acceleration. The model without re-acceleration is described by the conventional parameter 
set discussed previously with — 250 GeV and with a Merritt halo profile with a — 0.17 (the middle value in the 
acceptable range of values). These parameters are chosen because they generate a sufficiently large positron signal to 
explain the INTEGRAL excess _.34J. Fig. [5] shows our results for our canonical model with no re-acceleration. 

For a scenario with re-acceleration, we consider the same particle mass and profile and take va = 35 km/s. See 
Fig. [S]for the results of our canonical model with re-acceleration. An Alfven velocity of 35 km/s is reasonable, and 
we find that this value yields the smallest differences between the needed cross sections for the haze and HEAT, for 
a given (see Appendix [B|. The cross-sections differ at most by a factor of about 2 for va = 35 km/s (see Tables 
IVIIII and IIX [) . while for va = 20 km/s they more often than not differ by factors of 3 and higher (see Tables IVll 
and IIXP . Throughout the following discussion, the behavior of other models is very similar to that of the canonical 
models, unless otherwise noted. We refer the reader to the tables in Appendix IB] for a listing of the best-fit values for 
Po bl) the corresponding values for all models studied. 

The plots of the positron fraction include the HEAT data with error bars as calculated by the HEAT collaboration 
[l9| . The synchrotron radiation plots include the haze data as calculated by Dobler and Finkbeiner in [gH, based on 
WMAP measurements taken directly south of the Galactic Center. The data are in bins separated by 1 degree in 
latitude with a width of 20 degrees in longitude I G [—10, 10]. 

B. Positron Fraction 

With the results from PAMELA on the horizon, the positron production from dark matter is clearly an important 
avenue to pursue. However, the positron fraction contains numerous uncertainties, both in the background and in the 
shape of the DM signal, should there be one. In this section, we discuss a number of the different parameters and 
their effects on the positron fraction signal. 

1. Decay Channel Dependence 

The signals are noticeably different for energies close to for the two different decay channels, that is for the 
different mass ranges of the mediator (j). See Fig. [71 The bump in the positron fraction, characteristic of the hard 
XDM e^e~ spectra, has a larger maximum value for the direct decay channel. Moreover, in the muon channel the 
ratio starts its drop to the background value at a lower energy, so that the high energy fall-off is more gradual. These 
features can be understood by comparing the injection spectra for the two decay channels as shown in Fig. [2] The 
fiat distribution of the direct channel takes on a larger value at high energies than does the distribution of the muon 
channel. More e+e^ pairs produced at high energies give rise to a larger positron fraction at these energies. This also 
results in a smaller best-fit cross-section for the direct channel. 
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FIG. 5: Canonical models for va = km/s. 



2. Dependence 

As mentioned previously, the decay of the boosted particle directly to light fermions results in a hard positron 
spectrum. This gives rise to a characteristic enhancement in the positron fraction at energies above 5 GeV. The 
maximum in the positron fraction occurs at an energy that depends on the mass of x the decay channel; for 
example, for no re-acceleration and = 250, GeV the maximum occurs at ~ 140 GeV for the direct channel, and 
at ^ 50 GeV for the muon channel. The rapid, high energy fall-off of the positron fraction for the direct channel 
provides a way to determine . For the case with no re-acceleration, the high energy fall-off of the positron fraction 



occurs at ~ (see Fig. 8(a) I, while for the case with re-acceleration, the fall-off occurs at an energy larger than 



and that depends on the strength of the re-acceleration (see Fig. 8(c) ). Fig. [8] shows the dependence of the ratio on 
rriy- for both the direct and muon decay channels for both the no re-acceleration and the re-acceleration cases. 

As the mass of x increases, the best-fit cross-section correspondingly increases. This arises from the lower number 
density, = p^/m^, for fixed dark matter density. However, because higher energy particles can propagate farther, a 
given energy bin for HEAT is sensitive to a larger volume as one increases the energy range of the injection spectrum. 
This partially compensates for the decreased number density, but significantly larger cross sections are still required 
for higher DM masses. 
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FIG. 6: Canonical models for va = 35 km/s. 



3. Dependence on Halo Profile 



The positron fraction is fairly robust with respect to the halo profile. This is what we would expect, as the local 
positron fraction at high energies depends mainly on the e"'"e~ pairs produced nearby. Electrons and positrons lose 
energy very quickly through a variety of mechanisms as they travel through the Galaxy, so the fluxes of electrons and 
positrons that we measure locally at energies larger than 1 GeV are produced fairly nearby. Since the DM profiles are 
very similar at our location in the Galaxy, we expect the fluxes of DM-produced particles to be similar. Fig. [9] shows 
the positron fraction for all five halo models. The e~^e~ ratios are virtually indistinguishable, and this is true for all 
values of m^, though we emphasize that the cross-sections are different. As the Merritt parameter a increases, the 
cross-section increases. This is easily understood, as a is inversely related to the cuspiness of the profile. A smaller 
DM density results in fewer e'^e^ pairs produced, requiring a larger cross-section to make the fit to the HEAT data. 



4- Dependence on Alfven Velocity 

Fig. [TO] shows the positron fraction for the two different values of non-zero va studied, 20 km/s and 35 km/s. It 
appears that the re-acceleration parameter va has a strong affect on the positron fraction, and one might conclude 
that this is due to an underlying effect on the dark matter contribution to the ratio. However, the plots are misleading. 
As mentioned earlier, the effect of re-acceleration on the particle fluxes is small at energies above about 10 GeV, so 
the dark matter contribution to the positron fraction above 10 GeV is the same regardless of the value of va- The 
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FIG. 7: Positron fraction for direct decay (2me < < 2m^) and muon decay (2m^ < < 2m^o) channels. 



differences in tfie ratios at higfi energies as sfiown in Fig. [TO] are due to the differences in the bacfcgrounds at those 
energies (see Fig. |T]). The background positron ratio, which depends on the primary electron flux and secondary 
electron and positron fluxes, is determined by fitting the data at energies below about 7 GeV to the first five HEAT 
data points. Re-acceleration has a significant effect on particle fluxes at these energies, so to fit the data for different 
values of va, we must vary the primary electron flux through its power law index, the break rigidity of this index, and 
the normalization of the flux. (The secondary fluxes are also affected by re-acceleration, but we don't directly vary 
these through the inputs to GALPROP.) As the strength of re-acceleration increases, the power law index and flux 
normalization must increase to maintain the fit to data (see Table The changes in the background at low energies 
necessary to fit to the HEAT data come with corresponding changes in the background at high energies. It is the 
resulting differences in the background positron fraction at high energies that give rise to the differences in the total 
positron fraction at those energies. 



5. Implications for PAMELA 



Soon, results from PAMELA should provide data for high energy positrons up to ^ 270 GeV. It is somewhat natural 
that if the HEAT excess is, in fact, arising from dark matter, that PAMELA should see additional signals. Except 
for extremely light particles, the increased sensitivity and statistics of the PAMELA experiment should improve the 
understanding of this situation. Whether we can claim a signal as a discovery of dark matter we shall return to in a 
moment. However, even absent the HEAT signal, the haze motivates us to consider what might be seen at PAMELA. 
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FIG. 8: rrix dependence of positron fraction. 



Let us assume that the haze is, in fact, arising from dark matter annihilations as we describe here. A simple 
examination of the tables in Appendix [B] shows that the cross sections that explain the haze are in general comparable 
to or larger than the cross sections needed to generate a signal consistent with HEAT. Consequently, for a wide range 
of parameters, a strong signal would be expected at PAMELA, if the scenario as described explains the haze. 

We show representative data for PAMELA in Figs. [Tl] and [121 In Fig. [Tl] we consider two cases, = 50 and 
250 GeV, and present predictions for the PAMELA results after three years of data collection using the previously 
discussed backgrounds fit to the HEAT data. In Fig. [T^] we p resent similar predictions for = 250 GeV using 
different backgrounds fit to the preliminary PAMELA data For these backgrounds, the parameters defining the 
primary electron and proton injection spectra and the propagation parameters, e.g. the diffusion coefficient and the 
strength of re-acceleration, are listed in Table IIIII The backgrounds are in agreement with the local measurements of 
the electron and proton flux data as well as the B/C ratio data. 

It is clear that for the higher mass situation (see Figs. 



re-acceleration or not, a strong signal should be seen, 



11(c) 11(d) and [T2|) . regardless of whether one includes 



lich would be difficult to ascribe to uncertainties in the 



background. There appears to be some degeneracy in the data between lighter mass candidates and the decay 
channel considered. In particular, because the muon channel does not have a hard cutoff, it may be difflcult to 
establish the mass of the WIMP with any accuracy. Additionally, it may be difficult to distinguish electron from 
muon decay channels if the energy cutoff is above the PAMELA threshold. 

For lighter mass particles, the situation is not as clear. For instance, if we consider the situation without significant 
re-acceleration, it is impossible to claim discovery of a 50 GeV WIMP at PAMELA in this scenario, while with re- 
acceleration it is. This arises simply from the shape of the assumed primary electron spectrum, which gives rise to 
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FIG. 9: Positron fraction for all halo models for = 250 GeV. 
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FIG. 10: Va dependence of the positron fraction. 
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the background spectrum. For re-acceleration, the background is lower, so the dark matter component needed to fit 
the HEAT data is larger than that for no re-acceleration. Therefore, the signal in the re-acceleration case is larger. 
Discovery may still be possible for lighter mass particles in a scenario with no re-acceleration as measurements of 
other cosmic rays pin down additional astrophysical parameters, and as the low energy positron and electron spectra 
are established. 

Fig. [m shows our predicted PAMELA results for to^ = 250 GcV for two different values of m^, 212 MeV and 269 
MeV, the limits of the range of masses for which the decay into muons is the dominant decay channel. The decay 
though a lighter scalar mediator with the same value of Po 1"^!) results in a positron fraction with a larger peak 
value and a faster decline to the background value. These differences in signal are small for large values of m-^ and our 
backgrounds and cannot be distinguished by PAMELA. However, for smaller values of and a different background 
ratio these differences are more pronounced and possibly distinguishable. 

We want to emphasize that the background is dependent on the choice of propagation parameters, even within the 
constraints of local CR data. However, regardless of the differences in the background, we show that a significant 
signal can be seen at PAMELA for nix < 250 GeV with values of p§ {a \v\) of the the order 10~ 
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FIG. 11: Predictions for the PAMELA experiment results in the XDM scenario. 



C. DM Synchrotron Component 

1. Decay Channel Dependence 



While the different decay channels lead to quite distinct results for the positron fraction, the differences in the 
synchrotron radiation for the direct and muon channels are more subtle. See Fig. 1141 The muon channel produces 
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FIG. 12: Predictions for the PAMELA experiment results in the XDM scenario using the prehminaxy PAMELA data. 
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FIG. 13: Dependence of the predictions for the PAMELA experiment results on m^. 



TABLE III: GALDEF parameters for the background fits to preliminary PAMELA data 
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an intensity distribution that is more peaked in the center of the Galaxy than that produced by the direct channel. 
This leads to slightly better agreement with the haze data in almost all cases, particularly for larger values of m^. 
The best-fit cross-sections for the muon channel are larger than those for the direct channel for a specific and DM 
halo profile (refer eg. to Tables HVl and W} . Since the contribution to the number of high energy e+e~ pairs is larger 
for the direct channel annihilation process, one would expect more s ync hrotron radiation. The power spectrum of 
synchrotron radiation scales as the square of the boost of the particle [6^| , so it is consistent that the harder electron 
and positron spectra produce more radiation. 
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FIG. 14: 250 GeV XDM contribution to synchrotron radiation at 22.5 GHz for the direct and muon decay channels. 



2. Dependence 

Fig. [T5I shows the XDM contribution to the synchrotron radiation for various dark matter masses. The intensity 
distribution is more peaked in the center of the Galaxy for smaller values of for both production channels, leading 
to slightly better agreement with the haze data in almost all cases. This cuspiness arises simply because the lower 
energy particles do not propagate away from the Galactic Center as readily. As is the case for the positron fraction, 
the best- fit cross-section increases with increasing m^. Again, this can be understood as arising from the lower number 
density of dark matter particles, though this is somewhat compensated by the higher energy particles providing more 
energy for synchrotron radiation. 



3. Dependence on Halo Profile 



The shape of the intensity distribution of 22.5 GeV synchrotron radiation is relatively model- independent, with 
the notable exception of the isothermal model. See Fig. [161 The isothermal model produces a significantly flatter 
distribution. We investigated two additional halo profiles which we do not present here, the Evans profile and the 
alternative profile [131, both of which are included with the GALPROP code. The Evans model produced a DM 
synchrotron distribution that was even flatter than that produced by the isothermal model, while the alternative 
model produced a more peaked distribution that gave good agreement with the haze data for all m^. Ignoring the 
Isothermal model, we note that the N.F.W. profile requires a larger cross-section to fit the haze data. For example, 
the best-fit value of {a \v\) is 9.54 x lO-^^GeV^cm-^s-i for N.F.W. versus 6.72 x IQ-^eGeV^cm-^s-i for Merritt 
(a = .17) for the direct channel decay with — 250 GeV and va — 35 km/s. At small radial distances (r ^ 8 kpc), 
the N.F.W. profile gives a smaller value for the dark matter density than docs the Merritt profile (see Fig. [3]), so fewer 

e~ pairs are produced and thus less synchrotron radiation. 

Fig. 1161 shows the dependence of the XDM contribution to the synchrotron radiation on the Merritt parameter a. 
The intensity is more peaked in the center of the Galaxy for smaller values of a, and this results in better fits to the 
haze data. Additionally, the best-fit cross-section is smaller for smaller a. We understand this in terms of the effect of 
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FIG. 15: nix dependence of XDM contribution to synchrotron radiation at 22.5 GHz. 



the parameter on the shape of the Merritt profile; a profile with smaller a is more cuspy. The additional contribution 
to synchrotron radiation in the center of the galaxy for the cuspier profile reduces the cross-section necessary for a 
good fit to the haze. Note that the constant offset in our fitting (see Section llVEp compensates for the reduction at 
larger galactic angles. 



4- Dependence on Alfven Velocity 

Re-acceleration has a significant effect on the particle fluxes only below about 10 GeV, so the intensity distributions 
of synchrotron radiation for different values of va are identical. The effects of re-acceleration on the positron fraction at 
high energies came about from changes to the background ratio. Since the DM synchrotron component is independent 
of the background e'^e^ fluxes, the effects of re-acceleration are minimal. A comparison of the cross sections for 
different values of in Appendix [B] shows that the best-fit values are essentially identical. 



VI. CONCLUSIONS 



Dark matter remains an elusive component of the matter in the universe. Nonetheless, there are reasons for optimism 
that we may soon, if not already, have indications from astrophysics as to its nature. Hints from the WMAP haze 
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FIG. 16: XDM contribution to synchrotron radiation at 22.5 GHz for various halo profiles. 



and from excess positrons in the HEAT experiment both suggest new primary sources for high energy positrons and 
electrons. One explanation would be annihilating dark matter. 

We have considered the implications of "exciting" dark matter for these experiments and for the future results from 
PAMELA, and have found that the HEAT excess, as well as the haze, both arise naturally in this framework. Because 
the dark matter annihilates into a light state which then decays directly or indirectly to electrons and positrons, the 
resulting e'^e~ particles are highly boosted, and thus more relevant for experiment. In this way XDM falls into a 
broader framework of "secluded" dark matter [t^I where the annihilation products of the WIMP are not SM states 
themselves, but do couple to SM states. Some interesting avenues to consider include a broader range of states into 
which (p can decay and an expanded model framework. We will pursue these in future work. 

If local dark matter densities are high {po ^ 1 GeV cm~'^), the natural range of s-wave cross sections {{aann ~ 
2.5 X 10~^^ cm'^s"^) needed to explain the appropriate relic abundance also yields signals in agreement with the haze 
and HEAT, and with significant implications for PAMELA. It is important to note that in this scenario one needn't 
resort to large astrophysical boosts to explain HEAT, and that in most cases, the cross section which explains the 
haze has a significant signal at PAMELA. In particular, if one restricts oneself to the range of parameters which also 
explain the INTEGRAL signal (namely 400 GeV), if XDM explains the haze, a PAMELA signal is expected for 

the halo models we consider here. 

Many of the uncertainties associated with the signal arise from the uncertainty in the primary electron spectrum, 
i.e. its shape and normalization. Measuring the e+e" fluxes at PAMELA, as well as other cosmic ray observables, 
should reduce some of the astrophysical uncertainties, which should help to establish the significance of the signal 
beyond what we can discuss here. However, it is clear that PAMELA, with its measurements that are dominantly 
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sensitive to the local environment, will not be able to significantly distinguish between the halo models we considered 
here. In the event of a clear signal, however, PAMELA should be able to make some statement of the mass of the 
dark matter particle, even if just a lower bound. 

In conclusion, we have found that XDM provides a simple framework for explaining all existing astrophysical 
electronic anomalies: the low energy INTEGRAL 511 keV line from positron annihilation, the multi-GeV HEAT 
positron excess, and the possible multi-GeV electron and positron source generating the WMAP haze. Such signals 
all lie naturally within the same parameter range, and, in fact, the light mediator necessary for the success of XDM 
implies, in general, highly boosted annihilation products, which contribute to the significance of the haze and HEAT 
signals. Upcoming results from PAMELA should shed light on whether this scenario is realized in nature. 
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FIG. 17: Energy spectrum of electrons and positrons from muon decay. 

APPENDIX A: "MUON" DECAY CHANNEL SPECTRA 

The probability density for electrons and positrons created via the direct channel is a uniform distribution given by 
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The spectrum for the muon channel is a bit more involved. The distribution of injection energies for the final state 
electrons and positrons is the convolution of the spectrum of the positrons in the muon frame with the spectrum of 
the muons in the Galaxy, or DM, frame. The positron spectrum resulting from muon decay has been well-studied 
and is known as the Michel spectrum [Tlj . It is given by (ignoring radiative corrections) 
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where N is a normalization factor and Emax — the endpoint of the spectrum. It is shown in Fig. 1171 

The spectrum for the muons in the DM frame is identical to that for the positrons in the direct channel (see 

Eq. (jAip ) with the replacement rrie m^. The convolution of the two spectra was done numerically. The injection 

spectra for e~ pairs produced via both channels for = 0.1 GeV and = 0.25 GeV are shown in Fig. [2] 
With the following approximations, me — and 71 > 10 (see below), valid for the and x masses of interest to 

us, the injection spectrum for the muon channel can be calculated analytically. In this case the probability density is 

given by 
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APPENDIX B: BEST FIT PARAMETERS 
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TABLE IV: Table of bost-fit values of po (a \v\) and corresponding for the direct decay channel, — .1 GeV, for va ~ Q 
km/sec. We give separately the best-fit values for 

Po {<^Hz \v\) in units of lO'^'^GeV^ cm s to the haze data at 22.5 GHz and 
Po {o^HT \v\) to the HEAT data. We also give the corresponding X/fz = X^X^ for 29 data points of the haze and Xht = Y1X^ 
for 9 data points of the HEAT positron fraction. 
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500 
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110.42 
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8.77 



TABLE V: As table HV] but 



for the muon decay channel, = .25 GeV, for va = Q km/; 



sec. 



XlcLlO iViijLlfcJi 


Ivl^yKjG V ) 


'Z / 1 1 \ 
Po \(7HZ \V\} 
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Xhz 
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50 
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100 
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11.20 
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2.86 


7.37 
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43.60 
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11.23 


7.39 
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109.71 
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28.33 


7.40 



TABLE VI: As table ITVl but for the direct decay channel, 



= .1 GeV, for va = 20 km/ 



XlcLlO iViijLlfcJi 


Ivl^yKjG V ) 


'Z / 1 1 \ 
Po \(7HZ \V\} 
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500 
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7.36 



TABLE VII: As table |IVl but 



for the muon decay channel, = .25 GeV, for va = 



20 km/sec. 
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TABLE VIII: As table lIVl but for the direct decay channel, 



= .1 GeV, for va ~ 35 km/: 







'2 / \ \\ 
Po \0-HZ \V\} 

{ V 1 n~^^PoV^r'm~^G~^^ 
1 X iU vjrt; V Clll S ) 


Xhz 


Po {^HT \v\) 
I X lU Vjrc V Clli o ) 


Xht 


r^.t .\\ . 


50 


3.23 


683 


1.10 


8.81 




iUU 




DD / 


9 on 
z.yu 


7 17 
/ . 1 / 






1 9 nn 


71 fi 


lu.yo 


Q 1 1 
y. 1 1 




ouu 


'^U.oD 


7f^7 
/ / 


oD.Oo 


1 n 99 

lU. zz 




oUU 


yo. / o 


7sn 

/ oU 


Sfi ^i9 
oD.DZ 


1 n 

lU.Oo 


Isothermal 


50 


4.77 


799 


1.10 


8.50 




1 on 


D.Z ( 


091 
yz i 


9 s^n 

Z.oU 


o.y 1 






1 Q (\7 

iy.D / 


1 ni 7 

lUl ( 


1 n 7S 

lU. ( o 


Q n'^ 
y.uo 




ouu 


(^9 /in 


1 fVi^ 
lUoO 


O4.00 


1 n /in 

1U.4:U 




son 

oUU 


1 /1Q 1 /I 


1 n'i'i 

lUoO 


Q/1 


1 n Qi 
lu.y 1 


Mrt a = U.lo 




1 1 

i.ol 




U.o/ 


1 n /I T 




1 nn 

iUU 


9 "^n 

Z.oU 


oyy 


9 "^9 
Z.OZ 


S 1 1 
o. 1 1 




zou 


O.oo 


RQn 
oyu 


Q 1 /I 
y. 14 


8 7"^ 




ouu 


iD.4Z 


71 A 


•^n 1 1 
ou.ii 


Q 9n 
y.zu 




son 

oUU 


■^7 7n 
o / . / u 


I OD 


7n ^i.'x 


Q 9Q 

y.zy 


iVirt a = U.lY 


50 


2.21 


717 


1.01 


9.29 




1 nn 

iUU 


9 QS 


fi77 


9 1^/1 
Z.U4 


7/11 
/ .41 




zou 


o.yo 


71 Q 


1 n 9*^ 
lu. zo 


ft ss 
o.oo 




f^nn 
ouu 


97 fif^ 
Z ( .00 


7'^n 

( OU 


'X^ 7R 
OO. ( o 


Q 78 

y. / o 




snn 

oUU 




779 
/ ( z 


7Q 9fi 
/ y .ZO 


1 n V\A 

1U.U4 


Mrt a = 0.22 


50 


3.68 


680 


1.14 


8.63 




100 


4.86 


696 


2.90 


7.08 
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TABLE IX: As table ITVl but for the muon decay channel, — .25 GeV, for va ~ 35 km/sec. 
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